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ABSTRACT 

The crossing of the Galactic disk by a Globular Cluster could produce star for- 
mation due to gravitational focussing or compression of disk material. We report on 
simulations of the effect on disk material which reveal that the crossing can some- 
times cause local gravitational focussing of disk material. We also present the salient 
points of a little-known paper by Levy (2000), that shows that strong compression 
can result from the shock wave generated by GC disk crossing. The main thrust of 
our paper is a search for remnants of disk crossings by Globular Clusters. Using the 
gravitational potential of the Galaxy to locate the position of the most recent crossings 
of a subset of fifty-four Globular Clusters reveals that systematic errors and uncer- 
tainties in initial conditions limit the scope for unequivocal identification. From the 
subset of fifty-four, six possible search sites with the best constraints are retained for 
further scrutiny. Three of the six potentially promising search areas in the disk are 
from Globular Clusters NGC 3201, 6397 and NGC 6838, for which we cannot rule 
out some observed star associations observed nearby as being remnants. The three 
other of the six areas are too large to provide meaningful identification of remnants. 
Also, a possible remnant (open cluster NGC6231) is shown not to be due to Globular 
Cluster impact, contrary to a previous report. In a more wide-ranging screening of 
one hundred and fifty-five Globular Clusters we identify which Globular Clusters are 
compatible with being responsible for the formation of any of the Galaxy's five most 
prominent Star Super Clusters. 
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1 INTRODUCTION 

One of the less-studied interactions within our Galaxy con- 
cerns the effect of the passage of a Globular Cluster (GC) 
through the disk, on the disk itself. On average, one Glob- 
ular Cluster will cross the Galactic disk about every mil- 
lion years. While the Globular Cluster is in the vicinity of 
the disk, its gravitational field will attract stars and gas in 
the disk. The result could be a density increase sufficient to 
cause star formation. Whereas the impacts of High Veloc- 
ity Clouds with the disk have been examined in some detail 
(Comeron and Torra, 1994), as they cause gas compression 
and star formation (and are perhaps the origin of the Gould 
belt), the effect of GC impacts on the disk has received al- 
most no attention. The earliest mention of the effect of disk 
crossing by a Globular Cluster on the disk, is from Brosche 
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et al. (1991) who calculated the orbits of Globular Clusters 
NGC362 and NGC6218, and suggested observers examine 
their latest crossing points, to see whether any remnants 
may be present. Wallin, Higdon & Staveley-Smith (1996) 
developed the idea. Their main argument is that gravita- 
tional focussing of disk material, such as gas, towards the 
GC trajectory would occur. They concluded that large OB 
star associations would form, possibly some 30 Myr after 
the disk crossing, over a region of size ~10 pc. Their paper 
elicited little reaction, aside from two instances. The first 
was by Levy (2000), who examined compression of disk gas 
by the shock wave produced during GC travel through the 
disk. This paper appears never to have been cited, in spite 
of its importance to this question; we refer to it below. The 
second reaction to Wallin et al's paper was the suggestion by 
R. Rees and K. Cudworth, (reported in Wright, 2004) that 
the star association NGC 6231 in the Galactic disk is the 
result of such an interaction by Globular Cluster NGC6397, 
less than 5Myr ago. They reached this conclusion using po- 
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sitional and proper motion data, together with a code to 
calculate orbits in the Galactic potential. They found the 
point of impact of NGC6397 to correspond to the position 
of open cluster NGC6231, a finding we discuss later. 

The phenomena that take place in disk crossings by 
GCs are akin to those encountered in minor mergers and 
galaxy harassment, where a large galaxy interacts with a 
smaller galaxy at low and high relative velocities, respec- 
tively. Cox et al. (2006) predict that the global star forma- 
tion rate decreases with increasing mass ratio, so that merg- 
ers with mass ratio > 20 induce hardly any additional star 
formation in the more massive galaxy, above its quiescent 
level. In the case of GCs, the mass ratios are even larger, 
but relative velocities are lower and local effects in the disk 
may still be significant. Indeed, minor mergers have been 
proposed as the origin of nuclear starburst rings in some 
local galaxies (e.g. Mazzuca et al. 2006). 

The objective of this work is to examine whether, in 
principle, star formation, or at least disk material compres- 
sion could occur in a GC impact, and to look for possible 
remnants in the Galaxy. In section 2, we describe an ap- 
proach to modelling the interaction between a GC and the 
disk, which establishes how compression may occur via grav- 
itational focussing. This includes developing models of the 
Galaxy and GC potentials to calculate the trajectories of 
GCs and test masses in the disk. We then note an alternate 
mechanism for compression, namely shock wave generation. 
If gas is present, such a mechanism will be more effective at 
inducing star formation. In section 3, we apply the Galactic 
models to locate the most recent impact areas of GCs in the 
disk. Sections 4 and 5 use these location results to search 
for possible remnants in the Galactic plane. Finally, section 
6 provides conclusions. 

A preliminary version of this work was reported in 
Cropper & Vande Putte (2007). 



2 MODELLING 

We examine first the case of a gas-free disk, and consider 
the effect of a GC crossing, using numerical simulations of a 
series of test masses representing disk material encountering 
a Globular Cluster as it crosses the disk. These test masses 
are susceptible to the Galactic potential, and to the potential 
of the Globular Cluster, but not to one another. We later 
consider the case where gas is present at the impact site. 
First, however, we use the next two subsections to describe 
the methods we employ to calculate orbits in Galactic and 
GC potentials. 



2.1 The Galactic potentials 

We use four Galactic potentials for which straightforward 
analytical expressions are available. We denote these DA, 
FE, FL, and PA, as they appear in Dauphole & Colin (1995), 
Fellhauer et al. (2006), Flynn, Sommer-Larsen & Chris- 
tensen (1996), and Paczynski (1990), respectively. Others 
have used these expressions in studies of galactic dynamics, 
for example Kalirai et al. (2007), Law, Johnston & Majew- 
ski (2005), Pauli et al. (2003), and Dinescu, Girard & van 
Altena (1999). 



These potentials are all axisymmetric, so their expres- 
sions are simplest in cylindrical coordinates R, <j), z where R 
is the distance to the z axis, <j) is the azimuthal angle be- 
tween the x axis and the projection of the position vector 
on the xy plane, and z the distance above the x, y plane. 

The analytical characteristics of the potentials appear 
in Table 1. They consist of bulge, halo, and disk contribu- 
tions described by different models, as follows. 

The analytical expression for the potential from a Plum- 
mer (1911) sphere is: 

GM b 
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where r 2 — R 2 + z 2 . 

The Miyamoto-Nagai (1975) potential's analytical ex 
pression is: 
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The expression for the Hernquist (1990) potential is: 
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The logarithmic potential form is (Fellhauer et al. 
2006): 
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The logarithmic-arctan combination potential is 
(Paczynski, 1990): 
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Figure 1 shows a comparison between the potentials, 
radially along the disk, and as a function of z at R = 8kpc. 
The parameter values for the different cases appear in Table 
2. 

All calculations are undertaken for a fixed Cartesian 
Galactocentric right-handed coordinate system. The y-axis 
points from the Galactic centre to the Sun, and the z-axis 
points to the North Galactic Pole (NGP). The x-axis there- 
fore points in the direction of motion of the Local Standard 
of Rest (LSR). The initial position of an object results from 
its Heliocentric distance, and its celestial coordinates. The 
velocities of an object with respect to the LSR are deter- 
mined from its proper motion and radial velocity, and from 
the Solar motion (U Q = 10±0.4, V =5.2±O.6, W =7.2±O.4 
kms" 1 from Binney & Merrifield 1998, p628). In the LSR's 
coordinate system, the U axis points towards the Galactic 
centre, the V axis points in the direction of motion of the 
LSR, and the W axis points to the NGP. We transform ra- 
dial velocities and proper motions to motion in the LSR (U, 
V, W), using the formalism of Johnson & Soderblom (1987). 
The final transformation from velocities with respect to the 
LSR to velocities in a fixed Galactocentric system uses the 
LSR velocity given above for the four potentials, with an 
uncertainty of ±15 kms -1 suggested by Binney & Tremaine 
(1994, pl4). 

We prepared fortran95 computer codes to calculate 
orbits in the four potentials. We numerically integrate 
the equations of motion that arise from the potentials in 
cylindrical coordinates, using a fourth order Runge-Kutta 
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Table 1. Characteristics of selected potentials 



Potential 


Bulge 


Halo 


Disk 


Voire at 8kpc 
km s 


M total at lOOkpc 
10 11 M© 


Dauphole (DA) 


Plummer 


Plummer 


Miyamoto-Nagai 


225 


7.9 


Fellhaucr (FE) 


Hernquist 


log form 


Miyamoto-Nagai 


224 


9.3 


Flynn (FL) 


2xPlummer 


log form 


3xMiyamoto-Nagai 


222 


11.8 


Paczynski (PA) 


Miyamoto-Nagai 


log-atan form 


Miyamoto-Nagai 


220 


8.5 




Figure 1. Left panel, radial variation of potentials; right panel, vertical variation of potentials at R=8 kpc. Potential: DA, black; FE, 
red; FL, green; PA, blue. 



method. The codes require input of the initial position and 
velocity. We have verified the code by comparison with de- 
tailed numerical results from another institute (Turku Ob- 
servatory in Finland), and by comparisons with published 
orbits, and discussions with their authors. The codes gener- 
ate a series of files with the positions, velocities, total energy, 
^-component of angular momentum and elapsed time values. 
They also produce data on the position and timing of disk 
crossings, including the velocity at time of crossing. The user 
is able to specify the length of time and output times for the 
Runge-Kutta integration. All standard routines are from the 
NAG library 0. 

2.2 Addition of a moving potential 

The codes can also include a moving potential, such as that 
from a GC. In this case, the massive object experiences the 
potential of the Galaxy, whereas the test masses experience 
the potential of the Galaxy, as well as that of the massive 
object, but not of one another. We represent the GC's po- 
tential by a Plummer expression, where the parameter b in 
Eqn. (1) is derived from the mass of the Globular cluster 
and its central density po, using the expression from Binney 
& Tremaine (Eqn. 2-47b, 1994): 

1 www.nag.co.uk 



2.3 Gravitational focussing on a disk scale 

Wallin et al. (1996) make their case for star formation by 
considering phenomena at a disk scale of ~ O.lkpc, and we 
now investigate this in detail. 

Taking a set of fifty-four GCs described in more de- 
tail in section 3.1, we calculated the relative velocity, in the 
Galaxy, of each GC with respect to a cloud of disk mate- 
rial (test masses) on a circular orbit that meets the GC at 
disk crossing. We then model the cloud as 30,000 test masses 
placed at random in a box 0.1 kpc wide, 0.1 kpc deep, and 0.1 
kpc high, with the GC impacting the box along the normal, 
at the relevant velocity aimed at the centre of the box. We 
follow the test masses in the potential of the GC, until the 
GC leaves the box. Figure 2 shows the cloud of test masses 
in the box viewed along the GC trajectory for NGC5139 (u) 
Cen), the most massive Galactic GC, albeit here for only 
2,000 test masses, to avoid crowding the diagram. Gravi- 
tational focussing towards the centre is clearly visible, and 
is expected to be even larger when the gravity of the test 
masses themselves is accounted for. 

We repeated this experiment for a range of GC masses 
and velocities that cover the cases represented by our set of 
fifty-four GCs, and examined the compression as a function 
of GC mass/velocity, 50Myr after start of the experiment, 
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Table 2. Parameter values for the potentials 







Halo 




Disk — 






Bulge 








Dauphole 


M=1.3955xl0 10 
6=0.35 


M=6.9776xl0 n 
6=24 




M=7.9080xl0 10 
a=3.55 
6=0.25 




Fellhauer 


M=3.4xl0 10 
a=0.7 


1)0=186 km s — 1 
d=12 




M=10 n 
a=6.5 
6=0.26 




Flynn 


Af=3xl0 9 M=1.6xl0 10 
6=2.7 6=0.42 


v =220 km s _1 
d=8.5 


M=6.6xl0 10 
a=5.81 
6=0.3 


M=-2.9xl0 10 
a=17.43 
6=0.3 


M=3.3xl0 9 
a=34.86 
6=0.3 


Paczynski 


Af=1.12xl0 10 
a=0 
6=0.0277 


M=5xl0 10 

r h =e.o 




M=8.07xl0 10 
a=3.7 
6=0.2 





§Masses (M) are in units of Solar mass, halo velocities vq in kms -1 , and remaining parameters in kpc. The numbers of significant 
figures for parameter values reflect the situation in the original papers. The second disk mass in the FL potential is negative, but the 

disk and total densities are nevertheless positive everywhere. 
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Figure 2. View of box of 2,000 test masses after NGC5139 (w 
Cen) traverses it 



when all GCs have exited the box. Compression is shown in 
Figure 3, for three values of the parameter b in the Plummer 
potential (Eqn.l). Generally, compression increases with in- 
creasing mass, decreasing velocity, and is affected by the 
value of b. For small values of mass/ velocity, compression 
is higher for small b. As mass/velocity increases, the po- 
tential gradient is steeper, and the test masses are rapidly 
scattered away from the GC. For higher 6, the gradient is 
smaller (softer potential) and test masses spend a longer time 
in the vicinity of the GC track, leading to a ring of increased 
compression, due to an overlap of test masses moving to, and 
from the GC. 



2.4 Shock wave generation and material 
compression 

As modelled above, the test masses obey the collisionless 
Boltzmann equation, and we have shown that under such 
circumstances, gravitational focussing can occur. Even if the 
GC disk crossing does cause gravitational focussing, the per- 
turbation does not necessarily lead to star formation. This 
depends on the availability of molecular gas. Most of the 
Galaxy's molecular hydrogen is inhomogenously distributed 
and concentrated in a few thousand clouds of lOpc radius 
with masses >10 Mq (Binney & Tremaine, 1994), though 
sizes reach up to lOOpc with masses >6 xlO 6 Mq (Solomon 
et al., 1987). They are more prevalent in spiral arms, where 
they provide a fill factor of 0.1-0.5 (McKee & Ostriker, 
2007). 

In treating the impact on any gas in the disk, we need to 
include the additional collisional effects. The importance of 
this was first pointed out and investigated by Levy (2000). 
Levy treats the GC with an associated Plummer potential, 
impinging vertically on a gas layer. He uses a 2d hydrody- 
namical mesh code, and shows that a shock wave is set up 
by the GC. After dissipation of the shock wave, there re- 
mains a disturbed cylindrical region of diameter up to lkpc 
for a typical scale height of ~ lOOpc. Levy reports com- 
pression values in the shock wave of 2-30 for the different 
runs. If heating and cooling were allowed, these values could 
be larger. Levy also finds that compression increases with 
increasing mass, with decreasing velocity and with smaller 
b. We note that the compression values quoted here for the 
two mechanisms are broadly similar, with a range up to ~30. 
However, as the shock wave mechanism occurs in gas, it is 
expected to be more efficient at star formation. 



3 LOCATING IMPACTS WITHIN THE DISK 

Since both the collisionless and collisional treatments, and 
particularly the latter when gas is present, tend to increase 
the concentration of disk material, we now examine whether 
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Figure 3. Compression as a function of GC mass divided by GC 
velocity, for three values of b 



there is any observational evidence of star formation as a 
result of such impacts. 



3.1 The sample of Globular Clusters 

The most widely used Galactic GC catalogue is that of 
Harris (1996), which is regularly updated: the latest ver- 
sion, available at www.physics.mcmaster.ca/Globular.html 
is dated 2003. This version of the Harris catalogue lists one 
hundred and fifty-five GCs. The catalogue does not contain 
proper motion information, and such data are difficult to 
determine reliably. We therefore have narrowed our study 
to GCs with well-determined proper motion values from the 
subset of Allen, Moreno & Pichardo (2006 and 2008), ex- 
cept for NGC 6397, where we use more recent parameter 



values published since 2006 (Kalirai et al. 2007). The dis- 
tances from the Sun, radial velocities, and proper motions 
for the fifty-four GCs are given in Table 3. 



3.2 The location of impacts 

We used the four potentials and all fifty-four GC initial data 
to calculate their prior orbits. From this we have determined 
the location and time of the most recent crossing of the 
Galactic plane. Then we used the circular velocity at the 
GC Galactic radius for the potential in use to determine 
the present-day location of the impact points, knowing the 
time since the crossing. The resulting present-day positions 
of these crossings are shown in Figure 4. A few cross outside 
the visible Galactic disk, but most pass through it. Note that 
the four points that cluster around x — —90 kpc, y = 80 kpc 
in the left-hand panel correspond to Palomar 3. The right- 
hand panel concentrates on the Galactic disk. The average 
GC speed at time of crossing is ~280 kms -1 with a standard 
deviation of ~90 kms - . 



3.2.1 Effect of the potential model 

Systematic errors in impact location arise from our uncer- 
tain knowledge of the Galactic parameters, such as the ex- 
tent and mass of components making up the potential, and 
from more detailed configuration effects, such as the bar and 
spiral structure that are not reflected in the homogeneous 
bulges and disks. We neglect these latter effects as there 
are still uncertainties about the structure of the bar and 
arms. However we can gauge the effect of the former un- 
certainties in extent and mass of components by comparing 
all four axisymmetric Galactic potentials. For example, in 
Figure 4 (left-hand panel), the four points at x ~ 0, y >20 
kpc correspond to the crossing point of NGC 5024, as calcu- 
lated with the four potentials. There is a significant scatter 
of ~ 20kpc in the impact position. The crossing occurred 
~200Myr ago. Figure 5 shows a restricted part of the Galac- 
tic plane in the region of the Sun. This illustrates that the 
more concentrated impact points correspond to the more re- 
cent crossings, NGC 104 at ~47Myr, NGC 1851 at ~43Myr, 
NGC 3201 at ~5.4Myr, NGC 6397 at ~3.7Myr, NGC 6752 
at ~30Myr, and NGC 6838 at ~15Myr. 

We will concentrate our discussion on the six cases 
noted in Figure 5. They are selected because NGC 104/ 
1851/ 3201/ 6397/ 6752 impact points are the closest to us, 
and thus afford the best visibility. We add the impact point 
location of NGC 6838, as the path between it and us crosses 
a less dense region of the disk. In all six cases, the times 
since impact are short enough for massive stars still to be 
seen. Table 4 provides the times of crossing and the relevant 
distances from the Sun for these six GCs for the nominal 
values of initial conditions in Table 3. The uncertainties in 
the position of the impact points are discussed in Section 
3.2.2. The corresponding uncertainties in time since impact 
are shown in Table 4, and their origin are also discussed in 
Section 3.2.2. 
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Table 3. Proper motions, radial velocities, heliocentric distances, and masses of selected GCs 
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NGC6752 


-0.69±0.42 


-2.85±0.45 


3.8±0.19 


-32.1±1.5 


2.1xl0 5 


1, 12, 19 


NGC6779 


0.3±1 


1.4±0.1 


9.4±1 


-135.9±0.9 


1.5xl0 5 


1, 6, 12 


NGC6809 


-1.42±0.62 


-10.25±0.64 


5t 


174.9±0.4 


1.8xl0 5 


1, 12 


NGC6838 


-2.3=b0.8 


-5.1±0.8 


3.6±0.3 


-22.9±0.2 


3.0xl0 4 


1, 12, 20 


NGC6934 


1.2±1 


-5.1±1 


14.9±1.5 


-412.2±1.6 


1.6x10 s 


1, 6, 12 


NGC7006 


-0.96±0.35 


-1.14±0.4 


40±2 


-384±0.4 


2.0x10 s 


8, 12 


NGC7078 


-0.95±0.51 


-5.63±0.5 


9.5±1.1 


-106.6±0.6 


8.0x10 s 


1, 6, 12 


NGC7089 


5.9±0.86 


-4.95±0.86 


11.2±1.2 


-3.1±0.9 


6.9x10 s 


1, 6, 12 


NGC7099 


1.42±0.69 


-7.71±0.65 


7.3f 


-184.3±1 


1.6x10 s 


1, 12 


PAL3 


0.33±0.23 


0.3±0.31 


81.6f 


83.4±8.4 


3.3x10 s 


1 


PAL5 


-1.78±0.17 


-2.32±0.22 


9.9±2.2 


-55±16 


2.0xl0 4 


1, 7, 12 


PAL 12 


-1.2±0.3 


-4.21±0.29 


19.5=b0.9 


27.8±1.5 


l.lxlO 4 


9, 12 


PAL 13 


2.3±0.26 


0.27±0.25 


24.3±1.1 


-28±1.2 


5.4xl0 3 


10, 11, 12 



§Refs: 1, Dinescu et al., 1999; 2, McLaughlin et al., 2006; 3, Anderson & King, 2004; 4, Stetson, 1981; 5, Walker, 1998; 6, Dauphole & 
Colin, 1995; 7, Scholz ct al., 1998; 8, Dinescu et al, 2001; 9, Dinescu et al., 2000; 10, Siegel et al., 2001; 11, Cote et al, 2002; 12 Allen, 

Moreno, and Pichardo, 2006; 13, Alcaino ct al, 1997; 14, Szekely et al., 2007; 15, van de Ven et al, 2006; 16, Sarajcdini , 1992; 17 
Dinescu et al., 2003; 18 Kalirai et al., 2007; 19, Renzini et al, 1996; 20, Grundahl et al, 2002; 21, Allen et al., 2008; 22, Faulkner et al 
1991; 23, Cote et al., 1995; 24, Hartwick & Hesser, 1973; 25, Geislcr ct al., 1995; 26, Harris 1996; 27 Casetti-Dincscu ct al., 2007; 28, 

Alves ct al., 2001 
f No uncertainties on distance could be found for these clusters. 
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Figure 4. Present-day positions of latest crossing for the four potentials and fifty-four GCs in the Galactic plane. The Sun (red circle) 
is at x = kpc, y = 8 kpc. Potential: DA, filled black circles; FE, red stars; FL, green triangles; PA, blue open squares. The black circle 
of 15 kpc radius centered on the Galactic centre represents the extent of the visible disk (Sparke & Gallagher, 2000). The plot on the 
right shows an expanded region in the disk. 



Table 4. Present day distance D of last crossing point from Sun (kpc) and times T since latest crossing (Myr), for the four potentials, 
for the nominal kinematic parameters in Table 3. Uncertainties in position of impact point are provided in Figure 7, and discussed in 
Section 3.2.2. 







DA 




FE 




FL 




PA 




D (kpc) 


T (Myr) 


D (kpc) 


T (Myr) 


D (kpc) 


T (Myr) 


D (kpc) 


T (Myr) 


NGC104 


2.2 


48.8 ±2.1 


2.2 


46.7 ± 1.9 


2.3 


52.6 ± 1.9 


2.3 


49.1 ±2.3 


NGC1851 


2.5 


43.2 ±8.0 


2.1 


42.5 ±7.8 


2.5 


43.6 ±8.1 


2.5 


44.3 ±8.7 


NGC3201 


2.5 


5.5 ±0.4 


2.5 


5.4 ±0.4 


2.5 


5.5 ±0.4 


2.5 


5.5 ±0.4 


NGC6397 


2.6 


3.7 ±0.2 


2.6 


3.7 ±0.2 


2.6 


3.7 ±0.2 


2.6 


3.7 ±0.2 


NGC6752 


2.8 


29.4 ± 1.7 


2.8 


29.5 ± 1.6 


2.7 


35.3 ± 2.2 


2.8 


29.3 ±1.8 


NGC6838 


4.5 


15.6 ±0.5 


4.3 


14.6 ±0.5 


4.4 


18.6 ±0.7 


4.4 


14.7 ±0.5 



3.2.2 Effect of uncertainties in measured parameters 

Uncertainties (random errors) in the values of the Sun's mo- 
tion in the LSR, velocity of the LSR, distances, radial ve- 
locities, and proper motions affect the error in disk crossing 
location. We have therefore made a calculation that includes 
the effect of all these uncertainties. To appreciate the effect 
of these uncertainties on the impact sites of the six selected 
GCs, we assume that the values for Solar motion, LSR mo- 
tion, distances, radial velocities and proper motions are all 
normally distributed, with standard deviations equal to the 
quoted experimental errors. We then performed a Monte 
Carlo simulation for each of the six selected GCs, and 20,000 
trials per GC. Figure 6 shows the combined effect of the pa- 
rameter uncertainties and different potentials, at the 90% 
level. This also shows that the effect of the potential differ- 
ences can be as large as the effect of the orbital uncertainties, 
specifically in the case of NGC104. The combined results for 
the FE potential appear in Figure 7, as 90% and 68% prob- 
ability contour levels. We note that, as expected, the values 



of uncertainties in time since impact, and the size of the area 
encompassed by the 90% contours rank in the same order. 



4 SEARCH AREAS FOR REMNANTS OF 
IMPACTS 

Here we explore the six search areas identified above. 

4.1 Is NGC 6231 a remnant of an NGC6397 
impact? 

The only claim to have discovered an impact remnant - a 
cluster - in the disk relates to one of the six cases examined 
here: NGC 6397. Rees & Cudworth, cited in Wright (2004), 
suggested the site of the last impact of NGC 6397 to cor- 
respond to the open cluster NGC 6231 (/ = 343.46 deg, b 
= ±1.186 deg). The cluster has an angular diameter of 15' 
(Piatti, Claria & Bica, 1998). We calculate its impact time 
to be ~3.7Myr ago, consistent with Rees and Cudworth's 
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Figure 5. Restricted portion of Figure 4; unmarked positions 
relate to NGC 6752 



estimate of less than 5Myr. We can compare our time since 
last crossing at 3.7Myr with further estimates of the age of 
NGC 3201 in the literature. Piatti, Claria & Bica (1998) 
estimate it to be 5Myr. van den Ancker (2007) estimates 
the age to be ~3Myr. Sana et al. (2007) derive an age of 
at least lOMyr, based on examining the low-mass stars in 
the cluster, but they also indicate that low-mass star forma- 
tion was followed by a starburst 1-4 Myr ago. (Interestingly, 
van den Ancker (2007) concludes that the stellar population 
in NGC6231 is not incompatible with star formation being 
induced by GC impact on the disk, given that the cluster 
appears to have few low-mass stars, which he suggests may 
be the signature of a violent origin). Thus apart from Sana's 
determination, these ages are consistent with the last cross- 
ing time. However, we do not confirm the link established by 
Rees & Cudworth. Figure 8 shows the position of NGC 6231 
and the 90% level for the NGC 6397 impacts, at the time 
of impact 3.7 Myr ago. Also shown is the contemporaneous 
position of NGC 6231 (black square) and the points corre- 
sponding to that position ±3<r from that position, based on 
Sana et al.'s (2007) distance estimate of 1637±30 pc. Note 
that NGC 6231 is ~7pc in diameter. Clearly, NGC 6231 is 
not a remnant of the NGC 6397 impact, being ~0.6kpc away 
from the contour. Regrettably, the proper motions used by 
Rees & Cudworth were never published, preventing any fur- 
ther comparison. 



4.2 Impact points of NGC 3201, NGC 6397, and 
NGC 6838. 

The six areas and subtended fields of view from Earth en- 
compassed by the contours in Figure 7 are smallest for NGC 
3201, NGC 6397 and NGC 6838. As a result, these offer 
the best prospects for identifying remnants of crossings. For 
NGC 3201, the area is enclosed within a radius of 7 ° on the 
sky, and a distance range of ±22% from the nominal impact. 
For NGC 6397, these radii are 3 ° , with distance ±12% and 
for NGC 6838, they are 6 ° , and ±20%. 



8 - 



4 - 



- 




-8-4 4 

X (kpc) 

Figure 7. Present-day Galactic plane location of the 90 % (solid 
line), and 68 % (dashed line) contour levels for probability of 
impact, for NGC104(blue), NGC1851 (black), NGC3201 (pink), 
NGC6397 (red), NGC6752 (green), and NGC6838 (orange). Nom- 
inal impact positions are black squares. The Sun (red circle) is at 
x = kpc, y = 8 kpc, and only the FE potential is used here. 



We searched the CDS Simbad databas^EI centred on re- 
gions of sky corresponding to these areas, and restricted the 
search to object types "Star cluster, CI*", "Open Galactic 
cluster, OpC", and "Association of stars, As*" as defined in 
the database. This leaves 15, 33 and 42 star groups respec- 
tively, for NGC 3201, NGC 6387, and NGC 6838, and these 
are listed in Tables 5, 6 and 7. 

We searched the relevant literature for details of age 
and distance. Where this information is available, we noted 
for each star group in Tables 5, 6, and 7, the parameter that 
is incompatible, age (A) by a factor of at least two, and 
distance (D) outside the bounds given above. We conclude 
that of a total of fifteen star groups that lie in the direction 
of NGC3201's possible impact remnants, ten are definitely 
not remnants, on the grounds of incompatible age or dis- 
tance (Table 5). In the remaining five cases, it is possible 
that the groupings could result from the impact, pending 
further work on determining the properties (age, distance) 
of the associations. Examining the corresponding data for 
NGC6397 (Table 6), out of thirty-three candidates, eight 
can be excluded. The corresponding data for NGC6838 (Ta- 
ble 7) indicates that twenty five out of forty two candidates, 
can be excluded. This means that we are able to discount 
some GCs, and identify potential candidates, pending avail- 
ability of further data. 



2 http://cdsweb.u-strasbg.fr/, status March 2008 
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Figure 6. Present-day position of the 90% contour levels for the 6 GCs and 4 potentials (DA— black, FE— red, FL— green, PA— blue). The 
Sun (red circle) is at x=0 kpc, y=8 kpc. 
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Table 5. Star groups within the directions of NGC 3201 90% contours in Figure 7 



Identifier Type Basis for exclusion Refs 



C 0820-360 


CI* 


A, 


D 


1 


C 0815-369 


CI* 


A 




2 


C 0813-306 


CI* 


A, 


D 


1 


NGC 2546 


CI* 


A, 


1) 


1 


ESO 430-18 


CI* 


A 




1 


[KPR2005] 45 


CI* 


A 




3 


C 0812-362 


CI* 


A, 


D 


1 


C 0810-324 


OpC 


A, 


1) 


1 


C 0812-318 


OpC 


A, 


D 


4 


NGC 2588 


OpC 


A, 


D 


5 



Potential candidates are C 0805-322, [DBS2003] 16, C 0807-343, C 0809-318, DSH J0807. 1-3603, pending distance and age data 
Refs: 1, Kharchenko, 2005a; 2, Lindoff, 1968; 3, Kharchenko et al., 2005b; 4, Giorgi et al., 2007; 5, Baume et al., 2004. 



Table 6. Star groups within the directions of NGC 6397 90% contours in Figure 7 



Identifier 


Type 


Basis for exclusion 


Refs 


C 1715-382 


CI* 


D 


1 


[BDS2003] 97 


CI* 


D 


2 


[BDS2003] 98 


CI* 


D 


2 


[BDS2003] 164 


CI* 


D 


2 


[BDS2003] 99 


CI* 


D 


2 


C 1714-355 


CI* 


A, D 


4 


[DBS2003] 179 


CI* 


D 


5 


NGC6318 


OpC 


A 


3 



Potential candidates are [DBS2003] 120; [DBS2003] 121; [DBS2003] 122; [DBS2003] 123;MFSW VI; MFSW V; MFSW IV; C 1717 358; 
MFSW III; MFSW II; MFSW I; DSH J1715.7-3843; [BDB2003]G351.61+.17; [BDS2003] 119; [DBS2003] 118; C 1712-393; C 1727-370; 
CI VDBH 214; C 1720-378; C1721-389; C 1728-368; DSH J1713. 2-3942; C 1715-387; [DBS2003] 166; [DBS2003] 165, pending distance 

and age data 

Refs: 1, Piatti & Claria, 2002; 2, Bica et al., 2003; 3, Piatti et al., 2000; 4, Kharchenko et al., 2005b; 5, Borissova et al, 2005 



6.3 - 



5.9 - 



5.5 - 



5.1 - 




NGC 6231 



-1.7 -1.5 -1.3 -1.1 -0.9 -0.7 -0.5 

X (kpc) 

Figure 8. Position of open cluster NGC 6231 and of the 90% level 
of NGC 6397 impacts (FE potential), at the time of GC cross- 
ing. The triangle and circle represent the points corresponding to 
today's distance estimate ±3<r, respectively. 



5 COULD SUPER STAR CLUSTERS BE 
REMNANTS? 

An alternative approach to linking star formation to GC disk 
crossings is to work in reverse, looking for evidence of un- 
usual star formation, and asking whether this could be due 
to a GC disk crossing. For this we have concentrated on Su- 
per Star Clusters (SSCs). Portegies Zwart (2004) identifies 
five such SSCs: Arches, Quintuplet, NGC3603, Westerlund 
1, and Westerlund 2. Their distance from the Galactic centre 
and age appear in Table 8. 



5.1 Monte Carlo impact analysis 

We first look for progenitors amongst our sample of fifty- 
four GCs as we have already determined their crossings. It 
is possible to eliminate some of these, on the basis of the time 
elapsed since the last crossing. The time t since a crossing 
is given by t — tc + tsF + tA where tc is the time interval 
between disk crossing and the onset of star formation, tsF 
is the duration of star formation, and tA is the age of the 
cluster. The first time (tc ) is not well known, so three cases 
are considered here. Firstly, it is possible that this time is 
small compared to the age of the cluster, so that we can 
neglect it. Secondly, we adopt the value of 15Myr given by 
Lepine & Duvert (1994) for the case of star formation after 
a High Velocity Cloud impacts the disk. Finally, we take 
the value of 30Myr suggested by Wallin et al. (1996). For 
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Table 7. Star groups within the directions of NGC 6838 90% contours in Figure 7 



Identifier 


Type 


Basis for exclusion Refs 




CI* 


A, D 


1 


C 1 Q1 ^ -U1 Q4 


CI* 


A 


2 


[-Dl^/O^UUOJ -LOU 


CI* 


D 


3 




CI* 


D 


3 


[BDS20031 1 54 


CI* 


D 


3 


fBDS2003l 146 


CI* 


D 


3 


r"RFi i 39nn i il i ^9 

[DJJOZUUO] lOZ 


PI* 


D 


3 




CI* 


D 


3 


[rds20031 i 51 


CI* 


D 


3 


r"Rrt^9nn^l ia^; 

[_D_L7DZUUOJ ±^±0 


PI* 


D 


3 


W.D CllioLrei 111 VVOl 


PI* 


D 


3 


rRD t i9nn^i 14s 


CI* 


D 


3 


[RDS2003] 147 


CI* 


D 


3 


[RDS2003] 144 


CI* 


D 


3 


[BDS2003] 143 


CI* 


D 


3 


[BDS2003] 142 


CI* 


D 


3 


[BDS2003] 141 


CI* 


D 


3 


[BDS2003] 140 


CI* 


D 


3 


NGC6793 


CI* 


A, D 


1 


[BDS2003] 139 


CI* 


D 


3 


[BDS2003] 138 


CI* 


D 


3 


DSH J1930.2+1832 


OpC 


A D 


4 


DSH J1933.9+1831 


OpC 


A, D 


4 


NGC6802 


OpC 


A, D 


5 


C 1942+174 


OpC 


A, D 


1 



Potential candidates arc [BDS2003] 12; [BDS2003] 156; C 1926+173; DSH J1926.0+1945; [BDS2003]157; [BDS2003] 13; C 1926+147; 
DSH J1942. 3+1939; C 1916+156; DSH J1942.8+1530; C 1940+210; DSH J1912.0+1716; C 1922+136; DSH J1937.3+1841; CI Czcrnik 

40; DSH J1925. 2+1356; CI Alessi 57, pending distance and age data 
Refs: 1, Kharchenko ct al., 2005b; 2, Carraro et al., 2006; 3, Bica, et al., 2003; 4, Kromberg ct al., 2006; 5, Netopil et al., 2007. 



the star formation time, we adopt a value of tsF~2 xlO 5 
years, suggested by McKee and Tan (2002) , for massive stars 
(M > 8 Mo). Taking the value tc =30 Myr, and an average 
SSC age of tA =5 Myr, we are able to identify 24 GCs as 
not being SSC progenitors, because the time since impact is 
too long compared with 35Myr, using the FE potential. As a 
result, the list of fifty-four GCs to consider for Monte Carlo 
analysis is reduced to 30 candidates. A Monte Carlo analysis 
of these 30 candidates, similar to those above, using the FE 
potential, then shows further that the SSCs could not have 
been produced by any of these 30 GCs, with the exception 
of NGC362, 5897, 6093, 6316, 6522, 6584, and 6779, whose 
impact footprint is compatible with both the Arches and 
Quintuplet. Hence here too, some GCs can be discounted, 
while others may be potential progenitors. 



tential, using the relation V esc = ^J2.\<&(D)\ (Binney and 
Tremaine, 1994). If V> V eS c, then the GCs last disk crossing 
cannot be the event that triggered the formation of the SSC. 

We made the comparison for all GCs in the Harris cat- 
alogue, all five Super Clusters, the FE potential, and three 
times tc discussed above. Table 9 summarises the analysis. 
We can identify the GCs that are definitely not progenitors 
of the SSCs. As expected and evident in Table 9, the case 
with the smallest set of excluded GCs is that where tc =30 
Myr, and Table 10 lists the definitely excluded GCs. On the 
other hand, if star formation as a result of disk crossing is 
prompt ( tc ~0) then we can exclude all GC crossings as 
triggering the formation of both NGC3603 and Westerlund 
2. If this is generally the case, it may indicate that SSCs are 
not generated by GC disk crossings. 



5.2 Escape velocity analysis 

Having considered the fifty-four GCs, we then broaden our 
scope to all one hundred and fifty-five GCs in the Harris cat- 
alogue, and note that we can also screen for candidates by 
using a weaker method but which does not require knowl- 
edge of any proper motions or radial velocities. As the po- 
sition of the GC at the time of impact is the same as the 
position of the SSCs ancestral gas cloud, it is possible to cal- 
culate the minimum velocity V required for the GC to cover 
the distance from that position, to its present-day position 
within a time t . The escape velocity at the Galactocentric 
distances (D) in Table 8 can then be calculated from the po- 



6 CONCLUSIONS 

Globular Cluster interaction with a Galactic disk has been 
studied in detail from the viewpoint of the GC, see for ex- 
ample Gnedin, Lee & Ostriker (1999). On the other hand, 
the effect on the disk has received scant attention (Brosche 
et al. 1991, Wallin et al. 1996, Levy 2000). 

The suggestions in these papers is that GC disk cross- 
ings could cause star formation. We have therefore explored 
the process in more detail. We prescribe the Galactic poten- 
tial, and allow the GC to move within it. We also model the 
relative movement of a cloud of disk material on a circular 
orbit, as it encounters the GC crossing the disk, and show 
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Table 10. List of globular Clusters that are definitely not progenitors of the five SSCs for the case tc =30Myr. 



GC 


ssc 


GC 


SSC 


GC 


SSC 


AMI 


AON W1 W9 

rt, v*^, IN ■ VV J-, VV z 


IN ^jV^JUOt 


N Wl W9 

1 "I ■ V V -L > VV i 


IN V_I V _ | 


N Wl W9 

IN , VV -L , VV ^ 


AM4 


IV W1 W9 
IN , VV I, VV Z 


in vjv_^ouy^t 


A O NT Wl W9 

A, Vc^, IN, VV-L, VV Z 


Pall 


1ST W9 

IN , VV Z 


Arp2 


N, Wl, W2 


NGC5824 


N, Wl, W2 


Pal2 


A, Q, N, Wl, W2 


Eridanus 


A, Q, N, Wl, W2 


NGC6205 


W2 


Pal3 


A, Q, N, Wl, W2 


ESOSC06 


Wl 


NGC6229 


A, Q, N, Wl, W2 


Pal4 


A, Q, N, Wl, W2 


IC1257 


N, Wl, W2 


NGC6341 


N, W2 


Pal5 


N, Wl, W2 


IC4419 


Wl 


NGC6426 


N 


Palll 


N, W2 


NGC1261 


Wl 


NGC6715 


N, Wl, W2 


Pall2 


N, W2 


NGC1904 


W2 


NGC6779 


N, W2 


Pall3 


N, Wl, W2 


NGC2419 


A, Q, N, Wl, W2 


NGC6864 


N, W2 


Pall4 


A, Q, N, Wl, W2 


NGC4147 


N, Wl, W2 


NGC6934 


N, W2 


Pall5 


A, Q, N, Wl, W2 


NGC5024 


N, W2 


NGC6981 


N, W2 


Pyxis 


A, Q, N, Wl, W2 


NGC5053 


N, W2 


NGC7006 


A, Q, N, Wl, W2 


Rupl06 


Wl 


NGC5272 


W2 


NGC7078 


N, W2 


Terzan7 


N, Wl, W2 


NGC5466 


N, W2 


NGC7089 


N, W2 


Terzan8 


N, Wl, W2 




A=Archcs; ( 


3=Quintuplet; N= 


=NGC3603; Wl=Westerlundl; 


W2=Westerlund2 





Table 8. Distance from Galactic centre, and ages of Super Star 
Clusters 





Distance 
from Galactic 
centre 
D (kpc) 


Age 
Myr 


Rcfs. 


Arches 


~ 0.1 


2.5±0.5 


1 


Quintuplet 


~ 0.05 


4±1 


2 


NGC3603 


8±0.4 


~ 1 


3,4 


Westerlundl 


~ 9.4 


1-10 


5 


Westerlund2 


9.8±1.7 


~ 2.5 


6 



Refs: 1, Figer et al., 2002; 2, Figer, McLean, Morris, 1999; 3, 
Stolte et al., 2004; 4, Stolte et al., 2006; 5, Vrba ct al., 2000; 6, 
Rauw ct al., 2007. 



Table 9. Percentage of all Globular Clusters from the Harris 
catalogue that could have spawned the Super Clusters, based on 
their velocity being less than the escape velocity. 



t c = Myr t c = 15Myr t c = 30Myr 



Arches 20 76 93 

Quintuplet 34 80 93 

NGC 3603 16 75 

Westerlund 1 1 60 82 

Westerlund 2 17 73 



that focussing of disk material will occur for several of the 
GCs. The compression is on a scale of ~ lOpc, as suggested 
by Wallin et al. The extent of compression increases with 
GC mass and time spent by the GC in the cloud. Even if no 
gas is present at the time of interaction, the concentrated 
disk material can subsequently attract gas, leading possibly 
to star formation. 

Should gas be immediately available in the impact re- 
gion, a potentially more effective mechanism exists, pro- 
posed by Levy (2000), based on shock wave compression by 
the GC transiting through the disk. The star forming region 
in this case is estimated to be of kpc size. 

We have then examined whether there is any observa- 



tional evidence for star formation as a result of a GC disk 
crossing. The uncertainties in the parameters in our sub- 
set of fifty-four GCs, combined with the uncertainty in the 
potential model are such that the impact regions are gener- 
ally too extensive, subtending large angles on the sky. For 
some nearby GC crossings this is not the case, however, and 
for NGC 3201, NGC 6397, and NGC 6838, the subtended 
angles were smaller, as determined by Monte Carlo simula- 
tions. Using the CDS we identified a number of candidate 
star clusters and associations within the impact regions. We 
were also able to exclude a number of clusters and asso- 
ciations within the 90% contour on the grounds of age or 
distance. The only prior candidate for GC disk crossing star 
formation, the open cluster NGC 6231, lies well outside the 
90% contour of NGC 6397 impacts, in contradiction of the 
suggestion by Rees & Cudworth in Wright (2004). 

We then reversed the approach, and considered whether 
any of the Galaxy's five most prominent Star Super Clusters 
(Arches, Quintuplet, NGC 3603, Westerlund 1 and West- 
erlund 2) could have resulted from an impact by any of 
the sample of fifty-four GCs, using age arguments, and the 
Monte Carlo technique applied above. From this we found 
that only 7 GCs are candidates to provide impact site pro- 
genitors for the Arches and Quintuplet SSCs, with none for 
the other SSCs. As a final test, we considered all 155 GCs 
in the Harris catalogue and determined which SSCs could 
conceivably have resulted from a GC disk crossing. We used 
the criterion that the minimum speed to travel from a prior 
position of the SSCs to the current GC position should not 
exceed the Galactic escape velocity. This eliminated many 
of the larger sample of GCs, depending on the time elapsed 
between disk crossing and the onset of star formation. If the 
time between the disk crossing and the onset of star forma- 
tion is limited to a few Myr, it is possible to exclude any of 
the Galactic GC disk crossings as a potential cause of NGC 
3603 and Westerlund 2 formation. 
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